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Abstract 

We consider the decay of a massive particle under the complete or partial domination 
of the kinetic energy density generated by a quintessential exponential model and we 
impose a number of observational constraints, originating from nucleosynthesis, the 
present acceleration of the universe and the dark-energy-density parameter. We show 
that the presence of kination causes a prolonged period during which the temperature 
is frozen to a plateau value, much lower than the maximal temperature achieved during 
the process of reheating in the absence of kination. The decoupling of a cold dark matter 
particle during this period is analyzed, its relic density is calculated both numerically 
and semi-analytically and the results are compared with each other. Using plausible 
values (from the viewpoint of particle models) for the mass and the thermal averaged 
cross section times the velocity of the cold relic, we investigate scenaria of equilibrium or 
non-equilibrium production. In both cases, acceptable results for the cold dark matter 
abundance can be obtained, by constraining the initial energy density of the decaying 
particle, its decay width, its mass and the averaged number of the produced cold relics. 
The required plateau value of the temperature is, in most cases, lower than about 40 GeV. 
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1. Introduction 

A plethora of recent data [1, 2] indicates that the energy-density content of the universe 
is comprised of Cold (mainly [3]) Dark Matter (CDM) and Dark Energy (DE) with density 
parameters [1]: 

(a) J^cDM = 0.24 ±0.1 and (b) Ode = 0.73 ± 0.12, (1-1) 

respectively, at 95% confidence level (c.L). Several candidates and scenarios have been 
proposed so far for the explanation of these two unknown substances. 

As regards CDM, the most natural candidates [4] are the weekly interacting massive 
particles, x's- The most popular of these is the lightest supersymmetric (SUSY) particle 
(LSP) [5]. However, other candidates [6] arisen in the context of the extra dimensional 
theories should not be disregarded. In light of eq. (1.1), the x-relic density, ^^h^, is to 
satisfy the following range of values: 

(a) Om<n^h'^ and (b) n^h'^ <(}.13. (1.2) 

Obviously, the Q^h^ calculation crucially depends on the adopted assumptions. Ac- 
cording to the standard cosmological scenario (SC) [7], x's (i) are produced through thermal 
scatterings in the plasma, (ii) reach chemical equilibrium with plasma and (iii) decouple 
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from the cosmic fluid at a temperature Tp ~ (10 — 20) GeV during the radiation-dominated 
(RD) era (we do not consider in our analysis, modification to the Priedmann equations due 
to a low brane-tension as in ref. [8] and some other CDM candidates (see e.g. ref. [9]) 
which require a somehow different cosmological set-up). The assumptions above fix the 
form of the relevant Boltzmann equation, the required strength of the x interactions for 
equilibrium production (EP) and lead to an isentropic cosmological evolution during the 
X decoupling: The Hubble parameter is H (x with temperature T oc R^^ {R: the 
scale factor of the universe). In this context, the n^h"^ calculation depends only on two 
parameters: The x mass, and the thermal-averaged cross section of x times velocity, 
{(7v). If, in addition, a specific particle model is adopted, {av) can be derived from 
and the residual particle spectrum of the theory (sec e.g. refs. [6, 11]). Consequently, by 
imposing the CDM constraint - eq. (1.2) - some particle models (such as the Constrained 
Minimal SUSY Model (CMSSM) [10, 11]) can be severely restricted whereas some others 
(e.g. the models of refs. [12] and [13] which produce higgsino or wino LSP respectively) 
can be characterized as cosmologically uninteresting due to the very low obtained 

In a couple of recent papers [14, 15], we investigated model- independently (from the 
viewpoint of particle physics) how the calculation of Q^h^ is modified, when one or more of 
the assumptions of the SC are lifted (for similar explorations, see refs. [16, 17, 18, 19, 20]). 
Namely, in ref. [14] we assumed that x's (i') decouple during a decaying-massive-particle, 0, 
dominated era (and mainly before reheating) (ii') do or do not reach chemical equilibrium 
with the thermal bath (iii') are produced by thermal scatterings and directly by the (f) 
decay (which, naturally arises even without direct coupling [21]). The key point in our 
investigation is that the reheating process is not instantaneous [22]. During its realization, 
the maximal temperature, Tmax, is much larger than the so-called reheat temperature, Trh, 
which can be taken to be lower than Tp. Also, for T > Trh, -ff oc with T oc R~^/^ and 
an entropy production occurs (in contrast with the SC). As a consequence, in the context 
of this (let name it. Low Reheating) scenario (LRS) the n^^h'^ calculation depends also 
on Trhj the mass of the decaying particle, m^, and the averaged number of the produced 
x's, N^. We found [14] that, for fixed and {crv), comfortable satisfaction of eq. (1.2) 
can be achieved by constraining Trh (to values lower than 20 GeV), and N^. E.g., 
ri^/i^ decreases with respect to (w.r.t) its value in the SC [23] for low N^^^s and increases 
for larger A^^^'s [24, 25] (with fixed and (av)). Both EP and non-EP are possible for 
commonly obtainable {av)^s [26]. 

Another role that a scalar field could play when it does not couple to matter (in contrast 
with the former case) is this of quintessence [29]. This scalar field q (not to be confused 
with the deceleration parameter [3]) is supposed to roll down its potential undergoing 
three phases during its cosmological evolution: Initially its kinetic energy, which decreases 
as T^, dominates and gives rise to a possible novel period in the universal history termed 
"kination" [30]. Then, q freezes to a value close to the Planck scale and by now its potential 
energy, adjusted so that eq. (1.1b) is met, becomes dominant. In ref. [15], we focused on 
a range of the exponential potential [31] parameters, which can lead to a simultaneous 
satisfaction of several observational data (arising from nucleosynthesis, acceleration of the 
universe and the DE density parameter) in conjunction with the existence of an early 
totally kination-dominated (KD) era, for a reasonable region of initial conditions [32, 33]. 
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Table 1: Differences and similarities of ttie KRS witli ttie SC, LRS and tie QKS (ttie various symbols arc 
explained in sec. 2, the subscript I is referred to the onset of each scenario and est stands for "constant"). 

As a consequenee, during the KD era we have: H cc T'^ with T oc R^^. Under the 
assumption that the x-decouphng oceurs after the eommcncement of the totally KD phase 
- the assumptions (i) and (ii) were maintained - we found that in this scenario (let name 
it QKS) riy^h^ increases w.r.t its value in the SC [19] (with fixed and {av)) and we 
showed that this enhancement can be expressed as a function of the quintessential density 
parameter at the eve of nucleosynthesis, fiq^Tj^g). Moreover, values of Qq^z^^) close to its 
upper bound require {av) to be almost three orders of magnitude larger than this needed 
in the SC so as eq. (1.2) is fulfilled. It is obvious that the QKS although beneficial [34] for 
some particle-models [12, 13] can lead to an utter exclusion of some other simple, elegant 
and predictive particle models such as the CMSSM [19, 34, 15]. 

It would be certainly interesting to examine if the latter negative result could be evaded, 
invoking the coexistence of the two seenaria above (LRS and QKS) i.e. if a low TrH) which 
can assist us to the reduction of Q.^h? , can be compatible with a quintessential KD phase. 
To this aim we first investigate the dynamics of an oscillating field {(f)) under the complete 
or partial domination of the kinetic energy density of another field, q. We name this novel 
cosmological set-up KD Reheating (KRS). A similar situation has been just approximately 
explored in refs. [35, 36], under the name "curvaton reheating" for considerably higher 
scales (in their case q is restricted to drive quintessential [37] or steep [38] inflation and 
<j) is constrained so as it acts as curvaton [39]). The numerical integration of the relevant 
equations reveals that a prominent period of constant maximal temperature, TpL, arises 
(surprisingly, similar findings have been reported in ref. [40] for another cosmological set- 
up) . TpL turns out to be much lower than Tmax obtained in the LRS with the same initial 
energy density (the similarities and the differences between the various scenarios can easily 
emerge from table 1). On the other hand, the evolution of q is just slightly affected and 
so, it can successfully play the role of quintessence, similarly to the QKS. The resultant 
Q,^h'^ reaches the range of eq. (1.2) with (i") ~ (lO"'^ - 10"^) when TpL < Tp (type 
I non-EP), (ii") iV^ ~ when TpL ~ Tp (type II non-EP). On the other hand, EP which 
is activated for TpL > Tp requires a tuning of {av) in order that interesting ri^/i^'s are 
obtained. The required Tpl's are mostly lower than 40 GeV [41]. 

The framework of the KRS is described in sec. 2, while the analysis of the VL^h"^ 
calculation is displayed in sec. 3. Some numerical particlc-model-independent applications 
of our findings are realized in sec. 4. Finally, sec. 5 summarizes our results. Throughout the 
text and the formulas, brackets are used by applying disjunctive correspondence, natural 
units (^ = c = A;b = 1) are assumed, the subscript or superscript is referred to present-day 
values (except for the coefficient Vq) and In [log] stands for logarithm with basis e [10]. 
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2. Dynamics of Kination-Dominated Reheating 

According to the KRS, we consider the coexistence of two spatially homogeneous, 
scalar fields q and (/). The field q represents quintessence. The particle with mass 
can decay with a rate into radiation, producing an average number A'^^ of x's with mass 
m^, rapidly thermalized. We, also, let open the possibility that x's are produced through 
thermal scatterings in the bath. The system of equations which governs the cosmological 
evolution is presented in sec. 2.1 and a numerically robust form of this system is extracted 
in sec. 2.2. In sec. 2.3 we present the various observational restrictions and in sec. 2.4 we 
derive useful expressions, which fairly approximate our numerical results. 

2.1 Relevant Equations 

The cosmological evolution of q obeys the homogeneous Klein-Gordon equation: 

(a) q + 3Hq + Vg = 0, where (b) y = Foe"^*/""? (2.1) 

is the adopted potential for the g-field, , q [dot] stands for derivative w.r.t q [the cosmic 
time, t] and the Hubble parameter, H, is written as: 

^ " VSmp + Pq + P'i> + Pn with (b) = m^n^ and (c) pq = i^Q^ + y, (2-2) 



the energy densities of x and q correspondingly and = Mp/vStt (where Mp = 1.22 x 
10"*^^ GeV is the Planck scale). The energy density of radiation [(/>], p^ {p<j^-, and the number 
density of x, ?^x' satisfy the following equations [45, 24] (A<^ = (m^ — NyjUy^/m^): 

+ 3i?p0 + V^p^ = 0, (2.3a) 
p^ + 4Hp^ - T^p^ - 2m^{av) (^ - n'^f ) = 0, (2.3b) 
+ 3Hn^ + {av) (n^ - n^^"^) - F^N^p^/A^m^ = 0, (2.3c) 

where the equilibrium number density of x's, n^, obeys the Maxwell-Boltzmann statistics: 

n^{x) = j^:^rnl x^/^ e-V^P2(l/x), where x = T/m^, (2.4) 

g = 2 is the number of degrees of freedom of x and Pn{z) = 1 + (4n^ — l)/8z is obtained 
by asymptotically expanding the modified Bessel function of the second kind of order n. 
The temperature, T, and the entropy density, s, can be found using the relations: 

(a) Pn = ^9p*T^ and (b) s = —gs, (2.5) 

where gp*{T) [gs*(T)] is the energy [entropy] effective number of degrees of freedom at 
temperature T. Their numerical values are evaluated by using the tables included in 
micrOMEGAs [46], originated from the DarkSUSY package [47]. 

Finally, to keep contact with the LRS, we express the decay width of (p, F^, in terms 
of a temperature through the relation: 



^t> = 45 Mp - V 72 ^2-^) 

Note that the prefactor 5 is different from our choice in ref. [14] (4) with = Trh and 
several others (e.g. 6 [48] and 2 [18, 27]). Our final present choice is justified in sec. 2.4.2. 
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2.2 Numerical Integration 

The integration of the equations above can be reahzed successively in two steps: The 
first step concerns the completion of the KD reheating process, while the second one regards 
the residual running of q from the onset of RD era until today. 

2.2.1 First step of integration. The numerical integration of eqs. (2.1) and (2.3a)-(2.3c) is 
facilitated by absorbing the dilution terms. To this end, we find it convenient to define the 
following dimensionless variables [45, 18, 14] (recall that R is the scale factor): 

U = P4>R\ /r = PnR^ fi''^ = ni^'^R' and /, = qR:\ (2.7) 

Converting the time derivatives to derivatives w.r.t the logarithmic time [32, 33] (the value 
of Ri in this definition turns out to be numerically irrelevant): 

T=ln{R/Ri) =^ R' = R and R = Rie^ (2.8) 

eqs. (2.1) and (2.3a)-(2.3c) become (prime denotes derivation w.r.t f): 

Hf'^ = -V,R\ (2.9a) 

Hf'^ = -r^U, (2.9b) 

HR^fi^ = V^UR^ + 2m^{av) (/^ - ff) , (2.9c) 

HR^f^ = -{av) {fl - f'f) + T^NJ^R^A^m^, (2.9d) 

where H and T can be expressed correspondingly, in terms of the variables in eq. (2.7), as: 

1 , / -in \ 1/4 

VSR^m 



H = ^ JmJ^ + /|/2i?3 + VR^ + U + fn/R and T = ( " (2-10) 



The system of eqs. (2.9a)-(2.9d) can be solved from to ff 25, imposing the following 
initial conditions (recall that the subscript I is referred to quantities defined at f = 0) : 

UiO)Rf = p^^pI /r(0) = /^(O) = 0, g(0) = and 4(0) = ^/2^,{m^ / H^) . (2.11) 

Note that, unlike in the LRS [14], the numerical choice of /^(O) and g(0) is crucial for the 
result of We let the first one as a free parameter, while we determine the second 

one via m^, assuming that the (j) oscillations commence at Hi ~ [45, 36, 35]. Since 
we obtain the last condition in eq. (2.11). 

2.2.2 Second step of integration. When p^ S> p,j, the transition into the pure RD era has 
been terminated and the evolution of q can be continued by employing the formalism 
of ref. [15]. More precisely, we can define a transition point, t^, through the relation 
Pyi^^t) / P4>{^t) ~ and then, find the corresponding -r.p, via the formula [15]: 



^0 9p*_ f ^ 
\9s*) 



(a) pjf,) = p"^^f (^) e-4^T where (b) T=\n{R/R^). (2.12) 



The running of q can be realized from t^^ to 0, by following the procedure described in sec. 
2.1.3 of ref. [15] with initial conditions: 

q{T^) = q{T^) and q{T^) = {fJR^){i^). (2.13) 
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2.3 Imposed Requirements 

We briefly describe the various criteria that we impose on our model. 

2.3.1 KD "Constraint". We focus our attention on the range of parameters which ensure 
an absolute domination of the g-kinetic energy at f = 0. This can be achieved, when: 

(a) n\ = ng{o) = i with (b) ng = pg/{p^ + p^ + p^ + p^) (2.14) 

the quintessential energy density parameter. 

2.3.2 Nucleosynthesis (NS) Constraint. The presence of pq and p^ have not to spoil the 
successful predictions of Big Bang NS which commences at about T^^g = —22.5 correspond- 
ing to Tns = 1 MeV [49]. Taking into account the most up-to-date analysis of ref. [49], we 
adopt a rather conservative upper bound on Oq(T^g), less restrictive than that of ref. [50]. 
Namely, we require: 

(a) = flqiT^g) < 0.21 (95%c.l.) and (b) Trd > 1 MeV, (2.15) 

where Trd is defined as the largest temperature of the RD era and 0.21 corresponds to 
additional effective neutrinos species 6Ni, < 1.6 [49]. Wc do not consider extra contribution 
(potentially large [35]) in the left hand side (l.h.s) of eq. (2.15) due to the energy density of 
the gravitational waves [51] generated during a possible former transition from inflation to 
KD epoch [37]. The reason is that inflation could be driven by another field different to q 
and so, any additional constraint arisen from that period would be highly model dependent. 

2.3.3 Coincidence Constraint. The present value of pq, p^, must be compatible with the 
preferred range of eq. (1.1b). This can be achieved by adjusting the value of Vq. Since, this 
value does not affect crucially our results (especially on the CDM abundance), we decide 
to fix = pq/Pc to its central experimental value, demanding: 

(]0 = p0 = 0.73. (2.16) 

2.3.4 Acceleration Constraint. A successful quintessential scenario has to account for the 
present-day acceleration of the universe, i.e. [1], 

(a) - 1 < K;g(0) < -0.78 (95% c.l.) with (b) Wq = /2 - V)/{f /2 + V) (2.17) 

the barotropic index of the g-field. In our case, we do not succeed to avoid [33] the eternal 
acceleration {w^^ > — 1/3) which is disfavored by the string theory. 

2.3.5 Residual Constraints. In our scanning, we take into account the following less restric- 
tive but also not so rigorous bounds, which, however, do not affect crucially our results: 

(a) 10^ GeV <m^< lO^'' GeV and (b) < 1. (2.18) 

The lower bound of eq. (2.18a) is imposed so as the decay of (/) to a pair of x's with mass 
at most 500 GeV is kinematically allowed. The upper bound of eq. (2.18a) comes from 
the COBE constraints [52] on the spectrum of gravitational waves produced at the end of 
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inflation [37]. In particular, the later constraint impose an upper bound on Hi < 10^^ [15] 
which is translated to an upper bound on m^, due to our initial condition Hj = m^. Note 
that this constraint is roughly more restrictive that this which arises from the requirement 
for the thermalization of the 0-dccay products. Indeed, in order the decay products of the 
(j)-&eld are thermalized within a Hubble time, through 2^3 processes we have to demand 
^</> ^ 8 X 10^^ GeV [53]. The later is crucial so that eqs. (2.3a)-(2.3c) are applicable. 
Finally, the bound of eq. (2.18b) comes from the arguments of the appendix of ref. [24]. 

2.4 Semi-Analytical Approach 

We can obtain a comprehensive and rather accurate approach of the KRS dynamics, 
following the strategy of refs. [14, 15]. Despite the fact that our main interest is focused 
on the g-domination (p^^ < pq^) analyzed in sec. 2.4.3, we briefly review in sec. 2.4.2 
the dynamics of the decaying-^-domination {p^-^ > pq^) for completeness, clarity and better 
comparison. We first (see sec. 2.4.1) introduce a set of normalized quantities which simplify 
significantly the relevant formulas. 

2.4.1 Normalized quantities. In terms of the following dimensionless quantities: 

(a) -p^^^^=p^^Jpl (b) ^0 = 1^0^° and (c) q = q/V?>m,^. (2.19) 

eqs. (2.1) and (2.2) take the form (we use p° = 8.099 x IQ-^'^/i^ GeV^ with h = 0.72): 

(a) Q = Hq' and (b) HQ' + mQ + Vq = Q with (c) h'^ Pq + p^ + p^, (2.20) 

where (a) V = Vbe~^^^, H = H/Hq, Q = Q/ and (b) Pq = + V. (2.21) 

We do not present the normalized forms of the residual eqs. (2.9b)-(2.9d), since we do not 
use them in our analysis below. 

2.4.2 Decaying-<^Domination. When p^^ > p^j, the evolution of the universe at the epoch 
before the completion of reheating, T S> Trh, is dominated by p^ (transition to a q- 
dominated phase is not possible since pq decreases steeper than p^. Consequently, 

(a) H ~ p^l"^ and so, from eq. (2.9b) we obtain: (b) p^ = p^je"^*. (2.22) 

Substituting eqs. (2.22) in eq. (2.9c) and ignoring the last term in its right hand side (r.h.s), 
we can easily solve it, with result: 

Pn = 1^4>P^'^^ - ^'^^) with f ^ = T^/Ho. (2.23) 

The function p^{i) (in accordance with refs. [18, 14]) reaches at 

"^^ax - ln(1.48) = 0.39, a maximum value p^^^ ~ 0.14 f ^p^//^. (2.24) 

with corresponding T^ax derived through eq. (2.5a). The completion of the reheating is 
realized at f = fj^jj, such that: 

2 2 - 

Pr(^rh) = P<^(^rh) ^ ^RH - -3 In-f^P^j-^/^ (2.25) 

Equating the r.h.s of eqs. (2.22b) and (2.5a) for f = fj^jj and solving the resultant equation 
w.r.t we get eq. (2.6) with = Trh- We checked that the resulting prefactor 5 assists 
us to approach more satisfactorily than in ref. [14] the numerical solution of p^ = p^. 
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2.4.3 g-Domination. When > p^j)^, an intersection of pq with p,^ at t= fj^^ is possible, 
since pq decreases faster than p^. If this is realized, we obtain a KRS with partial q- 
domination (g-PD) whereas if it is not, we obtain a KRS with total (or complete) q- 
domination (g-TD). In either case, the KRS terminates at f = fj^^ which is defined as the 
commencement of the RD era (to avoid confusion with the case of sec. 2.4.2 we do not use 
the symbol fj^jj, any more). In particular, 

^RD = -^KR for 9-TD, or t^^ = t^^ for ^--PD, (2.26) 

where at ij^^ ["^^r] the transition from the KD [(?!)-dominated] to RD phase occurs. 

For q-TD and f < fj,^ or q-PD and f < fj^^, the cosmological evolution is dominated 
by Pq in eq. (2.21b), where the term V is negligible - see eq. (2.11). Therefore, 

(a) H ~ pq^/'^ and so, from eq. (2.20b) we obtain: (b) Pq = pq^e"^"^ . (2.27) 

Assuming for a while that eq. (2.22b) gives reliable results for the evolution even with 
g-domination (see sec. 2.4.3. b), we can achieve a first estimation for the value of fj^^, solving 
the equation: 

P<^(^K^) = ^ - ln(P9i/p0i) /3. (2.28) 

Obviously when < or, equivalently, Pqii^^) < P</'('^rd) '^^ obtain a KRS with 
g-PD. Let us assume that eq. (2.25) gives reliable results for fj^^ in the case of g-PD (see 
sec. 2.4.3.b). If we insert eqs. (2.28) and (2.25) in the first of the inequalities above or 
eqs. (2.27b) and (2.22b) in the second of the inequalities above and take into account that 
Hi = pqi^^"^ = rrifp/Ho, we can extract a condition which discriminates the q-TD from the 
q-PB: 

r^>2.5 p^Ho/m^ forg-TD, or < 2.5 p^Ho/m^, for g-PD. (2.29) 

Apart from the numerical prefactor, the condition above agrees with this of ref. [35]. Using 
the terminology of that reference, during the q-TD [q-FD] the (p-&eld decays before [after] 
it becomes the dominant component of the universe. In the following, we present the main 
features of the cosmological evolution in each case, separately. 

2. 4.3. a Total g-Domination. Substituting eqs. (2.27) into eq. (2.9b) we obtain: 

p^ ~ exp (-3t- ^c^^ e^*^ with c^^ = pq-^/^ = F^/m^. (2.30) 

The difference of the expression above from eq. (2.22b) is the presence of the second term in 
the exponent, which (although does not cause dramatic changes in the early p^^-evolution) 
is crucial for obtaining a reliable semi-analytical expression for the p^ evolution for f < f . 
Indeed, inserting eqs. (2.30) in eq. (2.9c) and ignoring the last term in its r.h.s, we end up 
with the following: 

'^R - Pch e"^* with Ir(t) = ^ dii exp ^4fi - ^c^^e^*'^ . (2.31) 
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For relatively small i (so as 4f 3> the integration of /r can be realized 

analytically and we can derive the simplified formula: 

Pn-\^.,P4-i (l ^^^PLf (2-32) 

where f^-^^ can be found by solving numerically the equation 4TpLf — lOc^^e^'^PLf /3. Prom 
eq. (2.32) we can easily induce that the function Pj^{t) takes rapidly (for f > fj^p^ = 1.5) a 
maximal plateau value: 

Prpl - ^c^0P<^i for 1.5<f<'fpLf, (2.33) 

Combining the previous expression with eq. (2.5a), we can estimate quite accurately the 
constant, plateau value of temperature: 

TPL ^ (7.5 c^^p^jTT^gp^iTpi^)) ^'^ . (2.34) 

When 4f < c^^e^'^/3, the first term in the exponent of eq. (2.31) can be neglected and 
the evolution can be approximated by the expression: 

~P^-\%iP<t>. e-4*[(e-4*PLf-l) 

+ \ (Ei(-^e3^) - Ei(-^e^*PLO) ] for ^ > (2-35) 

where Ei(2;) is the second exponential integral function, defined as Ei(2;) = — J^^e~^ /tdt. 
However, the contribution to of the terms in the second line of eq. (2.35) turns out to 
be numerically suppressed. Therefore, we can deduce that p^ for f > fpLf decreases as in 
the RD phase. 

Using eqs. (2.31), (2.30) and (2.27b) for solving numerically the equations: 

(a) Pr(v) = ^<^(v) (b) Pj,(Tkr) = P5(Tkr), (2.36) 

we can determine the points i^^ [^krI where p^ commences to dominates over p^ [p^] (see 
fig. 1 and table 2). In the present case, as we anticipated in eq. (2.26), the hierarchy is 
< f j^j^ and so, f can be identified as . 

The g-evolution can be easily derived, inserting eq. (2.27a) into eq. (2.20b) and ignoring 
the negligible third term in its l.h.s. Namely, 

qc^V2i (^ q =\/2) for t< t^^. (2.37) 

For T > Tjjp, we have H ~ p^/^. Inserting the latter into eq. (2.20b), we can similarly 
extract: 

q = +V2(l- e-^^-^KR)") for f_ < f (2.38) 



where q^^ = (?(fj^p). It is obvious from eq. (2.38) that q freezes at about -fj^p ~ Tj^^ + 6 to 
the following value: 

q^c^q^^ + V2 (^ q' = 0) for f^p < t. (2.39) 

Comparing these results with those derived in the context of the QKS (see ref. [15]), we 
conclude that the presence of ^ modifies just the value of Tj^p - due to the presence of 
eq. (2.31) - and so, it does not affect crucially the g-evolution. 
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2.4.3. b Partial g-Domination. When t< t^^, the [pq] evolution is given by eq. (2.27b) 
[eq. (2.30)]. Equating the r.h.s of these equations for f = i^^^, a more accurate result for 
fj^^ can be achieved than this obtained from eq. (2.28). However, the correction becomes 

more and more negligible as Pqi/P(f>i decreases (note that eq. (2.14) remains always valid). 

For T > Tj^^, the p^p [p^] evolution takes the form that it has in the decaying-(/>- 
doniinated era - see eq. (2.22b) [eq. (2.23)]. Namely, inserting eq. (2.22a) into eqs. (2.9b) 
and (2.9c) and integrating from fj^^ to f > fp,^, we arrive at the following results: 

P0 = P0(^kJ (2.40a) 
and -p^ = -p^{i^,) e-^(*-*K.) + h,-p,,'/\-^-'^ (l - e-^^^-^K.)/^) , (2.40b) 

where ptp{ij^^) [^^('^k,;,)] '^^^ derived from eq. (2.27b) [eq. (2.30)]. 

Employing the expressions above and eq. (2.27b) for solving numerically eq. (2.36b) 
[eq. (2.36a)], we can determine the points [^^b] where p^ commences to dominate over 
Pq [P4>] ^S- 2 and table 2). In the present case, as we anticipated in eq. (2.26), the 
hierarchy is i^^^ < f and so, i^^ can be identified as f . The resultant f^j^ approaches 
^RH obtained by eq. (2.25) as pq^/p^-^ decreases. 

As regards the (/-evolution, this obeys eq. (2.37) for i < i-^^. For t > fj^^, inserting 
eq. (2.22a) into eq. (2.20b) and ignoring the negligible third term in its l.h.s, we can extract: 

9 = 9"k, + (l - e-'(*-*K^)/') for T,, < f (2.41) 

where q^^ = q{ij^^) — A/21n (p^j/p^j) /3. It is obvious from eq. (2.41) that q freezes at 
about fj^p ~ T-^^ + 3 to the following value: 

q^ ~ q^^ + 2V2/3 (^ ^ = 0) for t^, < f. (2.42) 

Comparing these results with those derived in the previous case, we conclude that for g-PD, 
q takes its constant value during the ^-dominated epoch and that — q^^ is less than 
— ^j^H ™ q-TD. On the other hand, pq continues its evolution according to eq. (2.27b) 
until 3> Tj^p (note that ^'(fj^p) = but O('^kf) 0) where p^ reaches its constant value, 
Pqp = V{qp)- The point fp^ can be easily found (compare with ref. [15]): 

O'fe J/2 = Viq^) => fp, = Xqj2V3 - \n{Vo/pq^)/6. (2.43) 

3. Cold Dark Matter Abundance 

Our final aim is the Q,^h^ calculation, which is based on the well known formula [56]: 

^x = Px/'^c = (Wp[!)(%A)(^t)"ix ^ ^^x^' = 2.741 X 108 {f^/sR^){^)mJGeY, (3.1) 

where a background radiation temperature of Tq = 2.726 is taken for the computation 
of So and and ~ 25 is chosen large enough so as /^^ is stabilized to a constant value 
(with g's fixed to their values at T^). Based on the semi-analytical expressions of sec. 2.4, 
we can proceed to an approximate computation of f-^, which facilitates the understanding 
of the problem and gives, in most cases, accurate results. 
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We assume that x's are thermahzed with plasma (see sec. 2.3.5) and non-relativistic 
(m^ > T) around the 'critical' points f* or (see below). Our semi-analytical treatment 
relies on the reformulated Boltzmann eq. (2.9d). Note that, due to the prominent constant- 
temperature phase during the KRS, we are not able to absorb the dilution term of eq. (2.3c) 
by defining a variable Y = n^/s'^'' with i^s dependent on the form of T — R relation as, e.g., 
in refs. [56, 14, 15]. However, we checked that our present analysis based on eq. (2.9d) — 
see also ref. [18] - is quite generic and applicable to any other case. 

Since is stabilized to a constant value at large enough fj ^ 25, we find it convenient 
to split the semi-analytical integration of eq. (2.9d) into two distinct, successive regimes: 
One for < T < t^^ (see sec. 3.1) and one for fj^^ < f < f j (see sec. 3.2 [54]). In 
both regimes we single out two fundamental subcases: x's do (sec. 3.1.2 and 3.2.2) or do 
not maintain (sees. 3.1.1 and 3.2.1) chemical equilibrium with plasma. In the latter case 
for < f < f , two extra subcases can be distinguished: The type I and II non-EP. The 
conditions which discriminates the various possibilities are specified in sec. 3.1.1. Note that 
in our investigation we let open the possibility that x's remain in chemical equilibrium even 
after the onset of RD era. 

3.1 The Evolution Before the Onset of the RD Era 

During this regime, H can be sufficiently approximated by eq. (2.27a) [eq. (2.22a) and 
eq. (2.40a)] for q-TD or q-PD and f < i-^^ [for q-PD and f > fj^^]. Also, p^, involved 
in eq. (2.9d), can be found from eq. (2.30) [eq. (2.40a)] for q-TB or g-PD and i < i^^ 
[for g-PD and i > i^^]- Finally, T (involved in the computation of f^) can be found by 
plugging eq. (2.31) [eq. (2.40b)] for g-TD or q-PD and i < i^^ [for q-FD and i > t-^^] in 
eq. (2.5a). As regards the O^^/i^ calculation, we can distinguish the following cases: 

3.1.1 Non-Equilibrium Production. In this case, S> fx"^ for any f < fj^^ (type I non-EP) 
or fx <C fx^ for any t < f* (type II non-EP) . Let us consider each subcase separately: 

3. 1.1. a Type I. Obviously the realization of this situation (/^ S> fx'^ for any t < fj^^) 
requires Nx 7^ 0, (since if Nx = 0, the maximal possible value of fx is f^)- Such a 
suppression of fx^ can be caused if TpL ^ m^/20, as we can deduce from eqs. (2.7) and 
(2.4). Since - /^^^ ~ /2, eq. (2.9d) takes the form: 

HR^f^ = -{(jv) fl + V^NxP4>RV^4>m4>^ (3-2) 

which can be solved numerically from f = to f . In most cases (see fig. 4) the second 
term in the r.h.s of eq. (3.2) dominates over the first one and so, we can analytically derive: 

/r = 4V».)wi.h = (3.3) 

The integration above can be realized numerically. 

S.l.l.b Type II. When TpL ~ m-,j,/20, f^ is not strongly suppressed and so, the condition 
fx ^ fx^ can be achieved. Since — fx^"^ ~ ~fx'^i ^Q- (2.9d) takes the form: 



HR^f'^ = {av) + T^NxP^RVA^m^. 



(3.4) 
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Integrating the latter from f = to fj^^ we arrive at = /^^{i^j^) with: 

(a) /^^(T) = /^(f) + 4^(f), where (b) /^(f) = f di, {av) fx"^ / H (3.5) 

JO 

and is found from cq. (3.3). Wc observe that the integrand of /" reaches its maximum at 
f = f^, where the maximal ^-particles production takes place. Therefore, let us summarize 
the conditions which discriminates the EP from the non-EP of x's: 

!> 3 , type II non-EP (non-EPII) 

< 3 and (/^//^'') (f,) > 10, type I non-EP (non-EPI) (3.6) 
< 3 and if^/f^^) {t.) < 10, EP 

where the numerical values are just empirical, derived by comparing the results of the 
numerical solution of eqs. (2.9a)-(2.9d) with those obtained by the solution of eq. (3.2) or 
eq. (3.4). 

3.1.2 Equilibrium Production. In this case, we introduce the notion of the freeze-out 
temperature, Tp = r(fp) = x^m^ [7, 56], which assists us to study eq. (2.9d) in the two 
extreme regimes: 

• At very early times, when f <C fp , x's are very close to equilibrium. So, it is more 
convenient to rewrite eq. (2.9d) in terms of the variable A(t) = fxi^) ~ fx^{^) as follows: 

A' = -Z^-^' - (at;) A (A + 2f^^) /HR^ + T^N^p^R^ /HA^m^. (3.7) 

The freeze-out point Tp can be defined by 

A(Tp) = 6f f^^ir,) A(tp)(a(tp) + 2rii,)) = S^Sf + 2) f^'^^i,), (3.8) 

where 6f is a constant of order one, determined by comparing the exact numerical solu- 
tion of eq. (2.9d) with the approximate under consideration one. Inserting eqs. (3.8) into 
eq. (3.7), we obtain the following equation, which can be solved w.r.t fp iteratively: 



(ln/^^)'(i 



= -{av)5Fi5F + 2)/^'5(tp)/((5f + 1)HR' 
+ V^NxP^R^/{5f + l)iJ/^i(Tp)A<^m^ (3.9) 

->■(-/?)'(*) (-0) 

Normally, the correction to Tp due to the second term in the r.h.s of eq. (3.9) is negligible. 

• At late times, when f ^ Tp, ^ and so, — ~ Substituting this 
into eq. (2.9d), the value of at Tj^^p,, = fxi^Ko) '^^^ ^® found by solving eq. (3.2) 
from f = Tp until f = Tj^p, with initial condition fx{%) = ((^F + 1)/x''('^f)- However, when 

= or the first term in the r.h.s of eq. (3.2) dominates over the second, an analytical 
solution of eq. (3.2) can be easily derived. Namely = /^^(tj^d)' where: 

(a) 4^(T) = (/^(Tp)-i + JF(f))-' with (b) Mi)= Tdi,^. (3.11) 

J Tp 

The choice Sp = 1.0^ 0.2 provides the best agreement with the precise numerical solution 
of eq. (2.9d), without to cause dramatic instabilities. 
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3.2 The Evolution After the Onset of the RD Era 

During this regime, the cosmological evolution is assumed to be RD, and so, H = 
The evolution of and is sufficiently approximated by the following expressions: 

P0 = P0(^Ho) exp (^-3(f - - \ (^^y (e2(*-%o) - l) j (3.12a) 
and p^ = pjT^Je-^(*-*HD), (3.12b) 

where Trd corresponds to defined in eq. (2.26) and p^pii^^ is evaluated from cq. (2.30) 
[eq. (2.40a)] for g-TD [g-PD] while P^{t^^) is found from eq. (2.31) [eq. (2.40b)] for g-TD 
[^-PD]. In order to prove eq. (3.12a), we start from the exact solution of eq. (2.3a) [22, 7] 
which includes besides the terms of eq. (2.22b) an extra exponential term. We replace 
the involved temporal difference, in the latter term, by the corresponding temperature one 
(using the time-temperature relation [7] in the RD era) and by eq. (2.6), as follows: 



r</>(i - ^rd) - 




- l) (3.13) 



where in the last step we have used the entropy conservation law, eq. (2.8) and the fact 
that we do not expect change of Qs* between Trd and T^. Finally, T (involved in the 
computation of f^) can be found by plugging eq. (3.12b) in eq. (2.5a). 

As regards the ^^^h? calculation, we can distinguish the following cases: 

3.2.1 Non-Equilibrium Production. In this case, ^ for f > fj^^. This is the usual 
case we meet, when the evolution for f < Tj^^ has been classified in one of the cases 
elaborated in sees. 3.1.1 and 3.1.2. For f > Tj^q, the evolution obeys eq. (3.2) with 
P4> [Pr] given by eq. (3.12a) [eq. (3.12b)]. This equation can be solved numerically from Tj^^ 
until with initial condition f^i^m) ~ fx^ derived as we described in sec. 3.1. Under 
some circumstances an analytical solution can be, also, presented. Namely, 

/; = 4(*,) ''■'I' = (/x(%D)"'+^RD(i))"', where (3.14) 

(a) JMr, = £ ^ o. (b) J^if^ = £ « Sif. (3.15) 

Eq. (3.15a) is applicable for = or for (av) » [Tfj^N^/ Afj,'m(p){p^R^){Tj^j^) whereas 
eq. (3.15b) is valid when {av) f^^ < {T^NJ A^m^){p^R^)iT^^). 

3.2.2 Equilibrium Production. In this case, ~ f^^ for some f > Tj^^. This can be con- 
sidered as an exceptional case , since it can not be classified in any of the cases investigated 
in sec. 3.1 and is met for g-PD with very low Qq{Tj^g). The il^/i^ calculation is based on the 
procedure described in sec. 3.1.2. In particular, eq. (3.9) and (3.10) are applicable inserting 
into them eqs. (3.12a) and (3.12b). The limits of the integration of eq. (3.2) are from fp 
to fj in this case. Under the circumstances mentioned in sec. 3.1.2, an analytical solution 
can be obtained too. Namely, = /^(fj), with /^(t) given by eq. (3.11). 
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4. Numerical Applications 

Our numerical investigation depends on the parameters: 

A, m^, T^, N^, m^, (av). 

Recall that we use q{0) = and pq^ = (m^/ifo)^ (which is equivalent with Hi = m^jy) 
throughout. Also, Vq is adjusted so that eq. (2.16) is satisfied (we present the used Vq's in 
the explicit examples of figs. 1, 2 and 3). 

In order to reduce further the parameter space of our investigation, we make three extra 
simplifications. In particular, since A determines just the value of Wq during the attractor 
dominated phase of the g-evolution [15] and has no impact on the Q^^h? calculation we fix 
A = 0.5. Note that agreement with eq. (2.17a) entails < A < 1.15 [15]. Furthermore, 
since it is well known that, in any case, Q.^h? increases with (see, e.g., refs. [14, 15]) we 
decide to fix it also to a representative value. Keeping in mind that the most promising 
CDM particle is the LSP and the allowed by several experimental constraints and possibly 
detectable in the future experiments (see, e.g. ref. [55]) range of its mass is about (200 — 
500) GeV (see, e.g., fig. 23 of ref. [3]), we take = 350 GeV. 

Let us clarify once more that {av) can be derived from rriy^ and the residual (s)-particle 
spectrum, once a specific theory has been adopted. To keep our presentation as general 
as possible, we decide to treat m^^ and [av) as unrelated input parameters, following the 
strategy of refs. [14, 15]. We focus on the case {av) = a which emerges in the majority of 
the particle models (see, e.g., refs. [6, 12, 13] and [58]-[60]). We do not consider the form 
{av) = bx which is produced in the case of a bino LSP [10] without coannihilations [57, 59]. 
However, our numerical and semi-analytical procedure (see sees. 2.2, 3.1 and 3.2) is totally 
applicable in this case also with results rather similar to those obtained for {av) = a, as 
we showed in ref. [15]. 

The presentation of our results begins with a comparative description of the two types 
of g-domination (the g-TD and g-PD) in sec. 4.1 and of the various types of ^-production 
in sec. 4.2. In sec. 4.3, we investigate the behaviour of fl^h!^ as a function of the free 
parameters and in sec. 4.4 we compare the obtained O^/i^ in the KRS with the results of 
related scenaria. Finally, in sec. 4.5 we present areas compatible with eq. (1.2). 

4.1 Complete Versus Partial Domination of Kination 

The two kinds of g-domination, g-TD and g-PD, described in sec. 2.4.3 are explored 
in fig. 1 and fig. 2, respectively. Namely, in fig. 1 [fig. 2], we illustrate the cosmological 
evolution of the various quantities as a function of f for = 10^ GeV, = 30 GeV and 
logyo^j = 69.5 [logp^j = 77] (the inputs and some key outputs of our running are listed 
in the first [fourth] column of table 2) . We present solid lines [crosses] which are obtained 
by our numerical code described in sec. 2.2 [semi-analytical expressions as we explain in 
sees. 3.1 and 3.2] so as we can check the accuracy of the formulas derived in sec. 2.4. In 
particular, we design: 

• log Pi with i = q (black line and crosses), i = ^ (gray line and crosses), i = R (light 
gray line and crosses) versus f, in figs. l-(a) and 2-(a). In both cases we observe that Pq 
decreases more steeply than p^, and p^ remains predominantly constant. On the other 
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X = In R / R 

(c) 



x= InR/R 

(d) 



Figure 1: The evolution as a function of i for = IQ^ GeV, = 30 GeV and log p^j = 69.5 of 
the quantities: log pi with i — q (black line and crosses), i = 6 (gray line and crosses), i = R (light 
gray line and crosses) (a), T for pq^ = (m^t/Hof [pq — 0] (hold [thin] line and crosses) (b), log(sii^) for 
Pqi = (m^/Ho)^ [pq = 0] (bold [thin] Hue and crosses) (c), q and g' (d). The solid lines [crosses] are obtained 
by our numerical code [semi-analytical expressions]. 

hand, in fig. l-(a) [2-(a)], we observe that: (i) we obtain 2 [3] intersections of the various 
hnes, (ii) the hierarchy of the various intersection points is t^-^ < Tj^p [fj^^ < fj^p < t^r], (iii) 
at the point of the last intersection (fp^ = fj^p [^p^^ = T^p^]), we obtain {Pq/ P(j)){ijio) 3> 1 
[{Pq/P4>)i^RD) ^ 1] ^ expected from eq. (2.29). 

• T versus f, for pg-^ = (m^/i^o)^ [Pq = 0] (bold [thin] line and crosses) in figs. l-(b) 
and 2-(b) (obviously the thin lines correspond to a LRS with the same p^J. In both cases 
we observe that T rapidly takes its maximal plateau value, which is much lower than its 
maximal value obtained in the LRS - see eqs. (2.34) and (2.24). However, in fig. l-(b) the 
transition from the p^ < p^ to the > p^ phase, takes place at i^^ ~ 16.08 [fpjj ~ 12.06] 
(where a corner [kink] is observed on the bold [thin] line), whereas in fig. 2-(b), the same 
transition takes place practically at a common point t^^ 2± 17.47 for both the LRS and 
KRS where a slight kink is observed on both lines. This is expected since in fig. 2-(b) we 
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Figure 2: The same as in fig. 1 but with log = 77. 

obtain q-PD and so, for f > i^^^ = 14.3, the KRS and LRS give similar results. 

• log(si?^) versus t, for = (m^/i^o)^ [Pq = 0] (bold [thin] line and crosses) in 
figs. l-(c) and 2-(c). In both cases we observe that the initial entropy (si?^)(0) is much 
lower in the KRS than in the LRS. At the points where we observe a corner or a kink on 
the lines of fig. l-(b) [figs. 2-(b)], a plateau, which represents the transition to the isentropic 
expansion, appears in fig. l-(c) [figs. 2-(c)]. The appearance of the plateau observed on 
the bold and thin lines for q-TD - fig. l-(c) - occurs at different points (f^j^ ~ 16.08 and 
fj^jj 2± 12.06), whereas for q-PD - fig. 2-(c) - the same effect is realized at a common point 
since ~ 17.5 and f^^ ~ 17.82. This is expected, since for f > = 14.3 and g-PD, 
the KRS almost coincides to LRS (with the same p^p^). 

• q and q' versus f, in figs. l-(d) and 2-(d). In both cases we observe the period of 
the g-evolution according to eq. (2.37) with constant inclination (\/2) and the onset of the 
frozen field dominated phase (g' = 0). However, we observe that the frozen field phase 
commences much earlier in fig. 2-(d) (for g-PD) and q takes a value lower than the one in 
fig. l-(d) (for g-TD) - see eqs. (2.39) and (2.42). 
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4.2 Equilibrium Versus non-Equilibrium Production 



Figs. 


1, 3-(ai) 


3-(a2) 


3-(bi) 


2, 3-(b2) 


Input Parameters 


A = 0.5, m„ = l(f GeV. = 30 GoV. = 35 


G(V' 


log 


69.5 


73.82 


74.3 


77 


{av) (GeV-2) 


10-10 


10-10 


2 X 10-9 


1.8 X 10-9 




10-6 











Vb 


7 X 10" 


2.25 X 10^ 


1.18 X 10^ 


5.5 X 10"^ 


Output Parameters 


TpL (GeV) 


1.3 


16.05 


21.1 


100 




— 


— 


— 


14.3 


Tk4> (GeV) 








94.7 




16.08 


16.09 


16.1 


17.5 


T^R (GeV) 


1.3 


16.05 


21.1 


100 




22.1 


17.1 


16.6 


15.5 


Tkr (GeV) 


0.0035 


6.31 


13.9 


70.2 


PqJP<t>l 


1.3 X 10^^ 


6.4 X 10^1 


2 X 10^^ 


4.2 X 10^^ 


(Pg/P0)(^RD) 


2.5 X 10^2 


1.2 X 10^ 


11.84 


2 X 10-^ 


^('^^Ns) 


0.01 


3 X 10-9 


5.6 X 10"^° 


7 X 10-^5 




15.25 


15.7 


15.8 


15.3 




265 


22.7 


17.45 


4.6 







12.6 


0.0016 


5 X 10-9 








16.5 


18.5 








23.45 


23.86 




0.11 


0.11 


0.11 


0.11 


^x^^ sc 


1.87 


1.87 


0.11 


0.12 




1.78 


1.78 


0.1 


0.11 


^x^^lp<^=o 


979 


2.86 


0.12 


0.12 



Table 2: Input and output parameters for the four examples illustrated in 6gs. 3-(ai), (aa), (bi) and (ba) 
(see also figs. 1 and 2). 

The various kinds of x-production encountered in the KRS, are explored in fig. 3. In 
this, we check also the accuracy of our semi-analytical expressions (which describe the 
fy^ and evolution), displaying by bold solid lines [crosses] the results obtained by our 
numerical code (see sec. 2.2) [semi-analytical expressions (see sees. 3.2 and 3.1)]. The thin 
crosses are obtained by inserting T (given as we describe in sees. 3.2 and 3.1) into eq. (2.4). 

The inputs parameters and some key outputs for the four examples illustrated in fig. 3 
are listed in tabic 2. In particular, we present t's and the corresponding T's of the possible 
intersections between the various energy-densities. Comparing the relevant results, we 
observe that as p^-^ increases, Pqi/ptf,^ decreases (note that eq. (2.14) remains always valid), 
and so, {Pq/ P(j)){i^Yi) ^^('^ns) decrease too. At the same time, eventually approaches 
fjjj^ and becomes larger than this for logp,^; = 77, where g-PD is achieved (in the other 
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Figure 3: The evolution as a function of r of the quantities log{n^R'^) (bold line and crosses) and 
\og{n^^R^) + /„jq (thin line and crosses) for = lO" GeV, = 30 GeV and; (ai) logp^j = 69.5, = 
10"® and {av) = 10"^° GeV-^ (as) logp^j = 73.82, iV^ = and (av) = 10'^° GeV'^, (bi) logp^j = 
74.3, A^x = and (av) = 2 x 10"^ GeV"^ and (bs) logp^j = 77, A^x = and (av) = 1.8 x 10"® GeV"^ 
f/nf = 95 in tiie case (ai) and /„eq = eisewierej. The solid lines [crosses] are obtained by our numerical 
code [semi-analytical expressions]. In all cases, we extract O^'^^ =0.11. 

cases we have q-TD). In the same table we provide f*'s, derived from the maximahzation 
of the intergrand in eq. (3.5), and we apphed the criterion of eq. (3.6). In the case of EP, 
fp's and Xp's derived from eq. (3.9) are also given. In all cases, we extract Q^h^ = 0.11 
and we show the resultant ^^h^ in several other related scenaria (see also sec. 4.4). 

We design the evolution of log (bold line and crosses) and log + /„eq (thin line 
and crosses) as a function of i for = 10^ GeV, = 30 GeV and: 



log/ 



69.5, iV^ = 10"^ (/„cq = 95) and (av) = 10'^^ GeV"^ in fig. 3-(ai). In this 



case, we obtain q-TD (the evolution of the various energy densities is presented in fig. 1- 
(a)). The quantity turns out to be strongly suppressed due to very low TpL <C m^/20 
and so, we obtain S> for any f < 25. This is a typical example of non-EPI, where the 
presence of > is indispensable so as to obtain interesting Q.^h'^. Fixing = 10~^ and 
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Fig. 


Ranges of the Lower a;- Axis Parameters 


X-Pro- 

DUCTION 


N^ = 


= 10-^ 


= 10-6 


= 10-5 


4-(ai) 


lo.lo — i6.\)6 
73.94 - 74 


68.56 - 73.39 
73.4 - 74 


68.56 - 73.56 

73.57 - 74 


68.56 - 73.74 

73.75 - 74 


non-EPI 
non-EPII 


4-(a2) 


1 6.0 — 1 6.0 
73.61 - 74.5 


68.56 - 73.44 
73.45 - 74.5 


68.56 - 73.53 
73.54 - 74.5 


68.56 - 73.72 
73.73 - 74.5 


non-EPI 
non-EPII 


4-(bi) 


3.54 - 3.69 


i — O.i / 

3.18-3.28 
3.29 - 3.69 


(— u.y ) — 6 


^— U.y ) — u.byo 


non-EPI 
non-EPII 


4-(b2) 


3.4 - 3.54 


1 Q 1 7 

3.18-3.28 
3.29 - 3.54 


— u.y ) — o 


— u.yj — u.oyo 


non-EPI 
PP 

non-EPII 


4-(ci) 


3.88-4.1 


5-6.5 


5.8-6.8 


6.5 - 6.8 


non-EPI 
non-EPII 


4-(c2) 


3-4 
4-4.3 


3 - 4.2 
4.3 - 6.5 


3-4 
5.7-6.8 


3-3.7 
6.5 - 6.8 


EP 
non-EPI 
non-EPII 



Table 3: The type of x-pi'oduction for various ranges of the lower x-axis parameters and 's in fig. 4. 

adjusting p^j, we achieve fl^h'^ = 0.11. The bold crosses are derived by solving numerically 
eq. (3.2) for t < i^^^ (since {av) is rather low, eq. (3.3) is also valid) and from cq. (3.14) 
for < f < fj. The constant-/^ phase commences at f ~ 16.1, where the integrand of 
in eq. (3.3) reaches its maximum. 

• logp^j = 73.82, = (/„eq = 0) and {av) = 10-^° GeV-^ in fig. 3-(a2). In 

contrast with the previous case, we obtain a less efficient g-TD and so, TpL turns out to 
be significantly larger (TpL ~ m^/20). This is a typical example of non-EPII since at 
f* ~ 15.7, we get < f^. By adjusting p^j, we achieve ^}^h^ = 0.11. The bold crosses 
are extracted from eq. (3.5) for f < fj^^ and from eq. (3.14) for Tj^^ < t < fj. The onset of 
the constant-/;;^ phase occurs at f* 2± 15.7. 

• logp^j = 74.3, N^ = (/„eq = 0) and {av) = 2 x IQ-^ GeV-^ in fig. 3-(bi). In 
this case, we obtain a weak (since is very close to i^^^ , as shown in table 1) g'-TD with 
TpL > Tp. This is a typical example of EP before the onset of RD era, since for 15.8, 
we get > f^^ and < f^j-, = i^^. Jl^ft,^ = 0.11 is achieved by adjusting {(iv). The 
bold crosses are extracted by solving numerically eq. (2.9d) after substituting in it H and 
r as we describe in sees. 3.2 and 3.1. They could be, also, derived from eq. (3.11) for 
f p < f < fj^Q and from eq. (3.14) for fj^^ < i < if. 

• logp<^j = 77, Nx = (/„-<! = 0) and {av) = 1.8 x 10"^ GeV"^ in fig. 3-(b2). In this 
case, we obtain g-PD (the evolution of the various energy densities is presented in fig 2-(a)) 
and TpL 3> Tp. This is a typical example of EP after the onset of the RD era, since for 
f* ~ 15.3 we get > f^^, and fp ~ 18.5 > fj^p, = f^j^ ~ 17.5. Q^h^ = 0.11 is achieved by 
adjusting {av). The bold crosses are extracted similarly to the previous case. They could 
be also derived from eq. (3.14) for fp < f < fj . 
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Figure 4: Q^^ft^ versus log p^j (ai,a2), logT^ (bi, b2) and log (ci,C2) for fixed (indicated in the graphs) 
Tcj, and m^, p^j and m^, p^j and correspondingly, and various Nx's indicated on the curves. We take 
= 350 GeV and {av) = 10"" GeV"^ [{av) = 10"* GeV^] (ai,bi,ci [a2,b2,C2]). Tie solid lines 
[crosses] are obtained by our numerical code [semi-analytical expressions]. The CDM bounds of eq. (1.2) 
axe, also, depicted by the two thin lines. 
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4.3 ^x^'^ ^ Function of the Free Parameters 

Varying the free parameters, useful conclusions can be inferred for the behavior of 
Q^h^ ^^'^ the regions where each x-production mechanism can be activated. In addition, 
a final test of our semi-analytical approach can be presented by comparing its results for 
^^h^ with those obtained by solving numerically the problem. We focus on g-TD, since 
the results for g-PD are similar to those obtained in the LRS (see sees. 2.4.2 and ref. [14]). 

Our results are displayed in figs. 4 and 5. The solid lines are drawn from the results 
of the numerical integration of eqs. (2.9a)-(2.9d), whereas crosses are obtained by solving 
numerically eq. (2.9d) as we describe in sees. 3.1 and 3.2 (comments on the validity of 
eqs. (3.3), (3.5) and (3.11) are given, too). The running of f-^ after the onset of the RD 
era - see eqs. (3.11) and (3.14) - although crucial for the final result (especially for weak 
q-TD) does not alter the behaviour of the solution as a function of the free parameters. 
The type of x-production for the lower x-axis parameters and the various A^^^^'s used in 
fig. 4 is presented in table 3. From this and taking into account the obtained Tpl's in each 
case (see below), we can induce that non-EP (non-EPI [non-EPII] for N^- 7^ [iV^ ~ 0]) is 
dominant for TpL < m^/20, whereas EP is activated for TpL > 771^/ 20. 

In figs. 4-(ai, bi, ci) [figs. 4-(a2, b2, C2)], we take = 350 GeV and {av) = 10'^° GeV-^ 
[{av) = 10~^ GeV~^]. We design n^h'^ versus: 

• log/0(^j (or J7q(77j^g)) in fig. 4-(ai) and (32) for = 30 GeV, m<^ = 10^ GeV and several 
N^s indicated on the curves. We observe that: (i) fiq(Tj^g) increases as p^^ decreases (since 
Pqi/P4>i increases too) and so, a lower bound on p^j can be derived from eq. (2.15a) - 
note that TpL decreases with p^j (sec eq. (2.34)) and ranges between (0.8 and 24) GeV, 
(ii) Oi^/i^ decreases with p^j in fig. 4-(ai) {{crv) = 10~^'^ GeV~^) whereas it increases as 
p^j decreases (for large p^/s) and decreases with p^j (for low Pi^j's) in fig. 4-(a2) {{o'v) = 
10~^ GeV~^). The last observation can be explained as follows: can be mostly given 
by solving numerically eq. (3.2) see table 3. fy. increases with p^^ and so, when (av) is 
large enough {r-^ 10^^ GeV^^) the first term in the r.h.s of eq. (3.2) becomes comparable 
to the second one and the solution of eq. (3.2) can be exclusively realized numerically. For 
lower p^j's, eq. (3.3) can be used and so, decreases with p^. The latter behaviour is 
dominant for low (av) ~ 10"^° GeV~^ as in fig. 4-(ai). 

• logT,^ (or Qq{Tj^g)) in fig. 4-(bi) and (b2) for logp^j = 70, = 10^ GeV and several 
N^s indicated on the curves. We observe that: (i) Qq{Tj^^) increases with T<^ (or F^) since 
(j) decays more rapidly. Therefore, an upper bound on can be derived from eq. (2.15a) 

- note that TpL increases with (see eq. (2.34)) and ranges between (0.1 and 23) GeV, 
(ii) r2i^/i^ increases with (or F,^) and it turns out almost (<T?;)-independent for N^- 7^ 

- this is, because can be mostly given by eq. (3.3) as shown in table 3, (iii) Q^^'^ 
((jt;) -dependent and increases rapidly for N^- = - this is because can be extracted from 
eq. (3.5) as shown in table 3; decreases rapidly with (and Tpl) due to the exponential 
suppression of fy'^. 

• logm^ (or rig(Tj^g)) in fig. 4-(ci) and (02) for logm^ = 70, Tip = 300 GeV and 
several N^s indicated on the curves. We observe that: (i) Oq(Tj^g) increases with 
(since pq^ = (m<^/i?o)^ increases too) and so, an upper bound on m<^ can be derived 
from eq. (2.15a) - note that Tpl decreases as increases (see eq. (2.34)) and ranges 
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10" 10" 10'" 10' 10' 10' 10' ' 10" 10" 10" 10' 10' 10' 10' 



<ov> = a (GeV^) <ov> = a (GeV^) 

(a) (b) 

Figure 5: fl^h^ as a function of (av) = a for various N-^ 's, indicated on the curves, = 350 GeV, p^-^ = 
10™, = 10^ GeV and = 30 GeV [T^ = 3000 GeV] (a [b];. The sohd hues [crosses] are obtained by 
our numerical code [semi-analytical expressions]. The CDM bounds of eq. (1-2) are, also, depicted by the 
two thin lines. 

between (3.5 and 31) GcV, (ii) Q^h'^ decreases as increases for non-EP (see table 3) 
and increases with m,^ for EP (see fig. 4-(c2) and table 3). This is, because can be 
found from eq. (3.11a) for EP and it increases with or H (given by eq. (2.27a)) since Jp 
enters the denominator, whereas can be derived from eq. (3.3) [eq. (3.5)] for non-EPI 
[non-EPII] and it decreases when H increases. 

The dependence of fi^/i^ on (av) can be clearly deduced by fig. 5. We depict Jl^/i^ as 
a function of (av) = a for various AT^'s, indicated on the curves, = 350 GeV, = 
10'^°, = 10^ GeV and = 30 GeV [T^ = 3000 GeV] in fig. 5-(a) [fig. 5-(b)]. For these 
parameters we obtain q-TB with Tpl = 1.78 GeV [Tpl = 17.5 GeV] and nq{T^g) = 0.0017 
[ng{T^,) = 0.14] in fig. 5-(a) [fig. 5-(b)]. 

Obviously, due to very low Tpl, in fig. 5-(a) we obtain exclusively non-EPI. When the 
first term in the r.h.s of eq. (3.2) is comparable to the second one (usually for large {av)^s) 
Q^h? increases as (av) decreases, whereas when the second term dominates, O^/i^ becomes 
(crt;)-indcpcndcnt according to eq. (3.3). On the contrary, in fig. 5-(b), we obtain EP for 
(av) > 10-'^ GcV"^ and non-EPII for {av) < 10"^ GeV"^. We see that fi^/i^ increases 
as (av) decreases for EP, in accordance with eq. (3.11) - note that this is the well known 
behaviour in the SC. Also, O^/i^ decreases with {av) for non-EPII when the first term in 
the r.h.s of eq. (3.5a) is dominant (as shown in fig. 5-(b) for A^ = and for A^ = 10~^ and 
{av) in the range (10-9 - 10-^) GeV-2) whereas it remains ((7i;)-independent for non-EPII 
when the second term in the r.h.s of eq. (3.5a) is dominant (as shown in fig. 5-(b) for 
= IQ-'^ and {av) in the range (10"^^ _ ^q-S) GeV'^). 

Let us, finally, emphasize that the agreement between numerical and semi-analytical 
results is impressive in most of the cases. An exception is observed in fig. 4-(a2) for large 
p^j's where {Pg/P(j>){i^j^) < 50 and so, the adopted approximate formula for H in eq. (2.27) 
is not so accurate. The need for numerical solution of eq. (3.2) makes the discrepancy more 
evident than in fig. 4-(ai) where eq. (3.3) is everywhere applicable. 
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4.4 Comparison With the Results of Related Scenaria 

It would be interesting to compare Q^h'^ calculated in the KRS with that obtained in the 
QKS and the LRS, taking as a reference point the value obtained in the SC {pg = = 0), 
^x^^lsc- '^^^ relevant variations can be estimated, by defining the quantities: 



(a) An, 



^ and (b) An 



n 

^'X" I CD 



Isc 



XlCD 



sc 



n^h^l 



sc 



(4.1) 



where CD represents the condition which specifies the scenario under consideration: pg = 
for the LRS or = for the QKS. We restrict our analysis on the parameters used in fig. 4 
and we present the relevant results in tables 4 and 5. In table 4 we present il^/i^|sc and 
^^x\p(f,=o for several f2g(Tj^g)'s. For the same r2q(T^g)'s we arrange Af]^|p^=Q and AO^ in 
table 5. The corresponding values of TpL and Trh are also shown. Let us, initially, clarify 
the basic features of the n^^h^ calculation within the other scenaria. Namely, 

T, 
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in- 



dependent and so, it depends only on 
rriy. and {crv). Since these variables are 
fixed in figs. 4-(ai,bi,Ci) and figs. 4- 
(a2,b2,C2), we obtain two O^/i^lgf^'s 
presented in table 4. n^h^l^^ increases 
as (av) decreases (see also table 2). 



n^h^lp^^o is (P<^i, rn^ 



Table 4: ^l^h'^\sc s.nd ASl^^lp^^Q (for several Q,q{T^g)'s) 
for the parameters of figs. 4-(ai, bi, ci) and 4-{a2, b2, 02). 



independent and it exclusively depends 
on r2g(Tj^g) for fixed {av) and [15] 
(under the assumption that the onset 
of the KD phase occurs for T > m^). As a consequence, for several fixed r2g(77j^g)'s (see 
table 4), n^h'^\p^^Q takes a certain value for the figs. 4-(ai, bi, ci) and another for figs. 4- 
(a2,b2,C2). It is obvious that we obtain a sizable enhancement w.r.t n^h'^\g^, which in- 
creases with ilg(Tj^g) and as (av) decreases (see also the last line of table 2). 

• n^^h'^lp^^Q can be found by solving numerically [14] eq. (2.3a)-(2.3c) where H is given 
by eq. (2.2a) with Pg = 0. Note that coincides to Trh in the LRS. As we emphasized in 
ref. [14], the resultant n^hi^lp^^Q is -independent for Tmax > "m^ - see table 5, fig. 4-(ai) 
and (a2). So, in principle, Jl^/i^jp^^g is only (m^, T^, A'^)-dependent for fixed (av) and 
Since a variation of m,^ or changes ilq(Tj^g) for the KRS, we expect a change to 



m 



m 



X> "X''' \pg=0 



the corresponding n^h \pg^Q too. However, due to the fact that Trh 3> 
turns out to be m^-independent too - see table 5, fig. 4-(ci) and (C2). Finally, its N-^- 
dependence appears only for Trh < m^/iQ - see table 5, fig. 4-(bi) and (b2). We observe 
that n-^h'^\p^^Q turns out to be very close to n-^h'^\^^^ when Trh > m^/2{) (see also table 2) 
and mostly lower than n^h?'\^^, for Trh < m^/2Q - see table 5, fig. 4-(bi) and (b2). 

Comparing AJ7^ with AJl^jp^^Q we observe that contrary to the QKS (where n^h?\p^^Q 
exclusively increases with ng{T^^)) AO^ depends on the way that the Og(Tj^g) variation 
is generated. E.g., from table 5 we can deduce that An^ mostly decreases as ng{T^^) in- 
creases due to a decrease of p^j or an increase of and increases with ^^(Tj^g), when this 
is caused by an increase of T^. This can be understood by the fact that, in the two former 
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Table 5: AQ^ and TpL in the KRS and Afi^lp^^Q and Trh in the LRS for several Q,q{T^g)'s and the 
residual parameters of figs. 4-(ai, bi, ci) and 4-(a2, b2, C2). 

cases, TpL decreases whereas in the latter case, it increases. Obviously, the Oj^/i^-reduction 
with TpL is stronger when N-^ = 0. 

Comparing AO,^ with AJ7^|p^=Q we observe that: (i) TpL turns out to be much lower 
than TrH) (ii) increases with A^^ much more efficiently than Vt^h?\p^^Q, (iii) Ail^ can 

be positive in many cases (especially for Ny. 7^ and/or TpL ^ m^/20) in contrast with 
Ail^lp^^o which is mostly negative, except for large {av) and Trh = 20 GeV [25, 14], (iv) 
Q^/i^ approaches Q.yh?\p^^Q as Oq(T^g) decreases (for rig(T^g) < lO"^") - see table 2. 

4.5 Allowed Regions 

Requiring Vt^h? to be confined in the cosmologically allowed range of eq. (1.2), one 
can restrict the free parameters. The data is derived exclusively by the numerical program. 
Our results are presented in fig. 6. The allowed regions are constructed for N-^ = and 

= 10~^. In fig. 6-(ai,bi,Ci) [fig. 6-(a2, b2, C2)], we fixed = 350 GeV and {av) = 
10-^° GeV-2 [{av) = 10~^ GeV"^]. We display the allowed regions on the: 

• — log plane for = 10^ GeV, in fig. 6-(ai) and (a2). In the allowed regions of 
fig. 6-(ai) for A^ = [A^ = lO"*'] we obtain q-FB and EP for log p^^ > 74.8 [logp.^^ > 72.2] 
and g'-TD with non-EPII [non-EPI] elsewhere, while TpL ranges between (20 and 40) GeV 
[(0.85 and 1.5) GeV]. Since increases with p^^ - see eq. (3.3) [eq. (3.5)] for non-EPI 
[non-EPII] and eq. (3.11) for EP - the upper [lower] boundaries of the allowed regions 
come from eq. (1.2b) [eq. (1.2a)]. The lower right limit of the allowed regions comes from 
eq. (2.15a). As increases Oj^/i^ approaches n^h'^\p =0 ^^'^ it becomes equal to 0.09 at 
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Figure 6: Regions allowed byeq. (1-2) on theT^— logp^j pJane (ai, 82) form^ = 10 GeV, logm,^— logp,^j 
plane (bi, b2) forT^ = 30 GeV and = 300 GcV, and T^ — \ogm^ plane (ci, C2) for logp^j = 70. Wc take 
N^ = OorN^ = 10-^ = 350 GeV and (av) = 10"'" GeV"^ [{av) = 10'^ GeV"^] (ai, bi, ci [a2, b2, C2]). 
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the upper left bounds of the allowed regions (possible further reduction of reduces also 
Qy^h"^ which turns out to be independent, any more). In the dark grey [grey] allowed 
regions of fig. 6-(a2), we obtain q-TD with EP [non-EPI] while TpL ranges between (18 and 
400) GeV [(1 and 15) GcV]. Due to the increase of the released entropy which is caused by 
the increase of p^^, ^^^^ decreases as p^j increases and so, the upper [lower] boundary of 
the dark grey and the almost horizontal part of the grey allowed region come from eq. (1.2a) 
[eq. (1.2b)] (note that eq. (3.3) is not applicable in this part of the grey allowed region). 
On the contrary, eq. (3.3) can be applied for the almost vertical part of the grey allowed 
region and so, its left [right] boundary comes from eq. (1.2a) [eq. (1.2b)]. The lower right 
limit of this region is found from eq. (2.15a) whereas the upper one is just conventional. 

• logm^ — logp^^jj plane for = 30 GcV and = 300 GcV, in fig. 6-(bi) and (b2). In 
the allowed areas we obtain gr-TD. In the dark and light grey areas of fig. 6-(bi) [fig. 6-(b2)] 
we have non-EPII [EP]. The dark grey regions are A'^-independent since the first term in 
the r.h.s of eq. (3.5a) [eq. (3.2)] dominates over the second one for non-EPII [EP]. In fig. 6- 
(bi) the upper [lower] boundary of the allowed regions comes from eq. (1.2b) [eq. (1.2a)], 
whereas in fig. 6-(b2), the origin of the boundaries is the inverse. This is because higher /J^j 's 
result to hither Tpl's and so, Trd's. This has as a consequence that the non-relativistic 
reduction of turns out to be less efficient, thereby increasing ^^h^ for non-EPII. For 
EP, the same effect causes mainly an increase to the released entropy and so, a reduction 
to U^h?. Note, also, that since c^^p^^ remains constant along the boundaries of the dark 
grey areas TpL remains also constant - see eq. (2.34) - and is equal to 16 GeV [18 GeV] 
in fig. 6-(bi) for = 30 GeV [T^ = 300 GcV] and (32 - 38) GeV in fig. 6-(b2). The 
upper [lower] limits of the dark grey or grey areas (in the upper right [lower left] corners) 
of these figures correspond to the upper [lower] bound of eq. (2.18). In the grey areas of 
fig. 6-(bi) and (b2), we obtain non-EPI. Eq. (3.3) is apphcable in the almost vertical parts 
of these areas and so, the left [right] boundary of the allowed regions comes from eq. (1.2b) 
[eq. (1.2a)] (TpL ranges between (3 and 17-18) GcV). Eq. (3.3) is not apphcable in the left 
upper branch of the area in fig. 6-(b2) where the origin of the boundaries is the inverse 
(TpL = (1 — 7) GcV). The lower bound of the almost vertical part of these areas come 
from eq. (2.15a). Note, finally, that for = 30 GeV and {av) = 10"^ GeV'^ - fig. 6-(b2) 
- we are not able to construct allowed area for iV^ = 0. This is because TpL ^ m^/20 
and possible increase of p^j (which could increase Tpl) leads to Q,^h? close to Vt-^h?\p^^Q 
which is lower than 0.09 due to large {(tv). However for N-^ = 10~^, Oj^/i^ increases to an 
acceptable level. 

• Tff, — logm<j plane for logP(^j = 70, in fig. 6-(ci) and (02). In the allowed regions of 
these figures, we obtain (/-TD, besides the part of the grey area for log(m<^/GeV) < 5.35 
where we get g-PD. For 7^ we have non-EPI, whereas for = 0, we take non-EPII, 
except for the lower part of the light grey area (for \og{m^/GeV) < 3.3) where the EP is 
activated. As induced from eq. (3.3) [eq. (3.5)] for non-EPI [non-EPII], decreases as 

(and so, p^J decreases. Therefore, the upper [lower] boundary of the allowed regions 
comes from eq. (1.2a) [eq. (1.2b)]. The upper bound on the right corners of the allowed 
regions is derived from eq. (2.15a). The lower limits of the light grey areas and these of 
the lower left corner of the grey areas are extracted from the lower bound of eq. (2.18a). 
In the grey [light grey] areas Tpl ranges between (0.1 and 5) GeV [(16 and 20) GeV]. 
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5. Conclusions 

We studied the decoupling of a CDM candidate, x, in the context of a novel cosmo- 
logical scenario termed KRS ("KD Reheating"). According to this, a scalar field (j) decays, 
reheating the universe, under the total or partial domination of the kinetic energy density 
of another scalar field, g, which rolls down its exponential potential, ensuring an early KD 
epoch and acting as quintessence today. We solved the problem (i) numerically, integrating 
the relevant system of the differential equations (ii) scmi-analytically, producing approxi- 
mate relations for the cosmological evolution before and after the onset of the RD era and 
solving the properly re-formulated Boltzmann equation which governs the evolution of the 
X-number density. Although we did not succeed to achieve general analytical solutions in 
all cases, we consider as a significant development the derivation of a result for our problem 
by solving numerically just one equation, instead of the whole system above. 

The model parameters were confined so as Hi = m<^ and $7^ = 1. The current ob- 
servational data originating from nucleosynthesis, acceleration of the universe and the 
DE density parameter were also taken into account. We considered two cases depending 
whether (j) decays before (g-TD) or after ((/-PD) it becomes the dominant component of 
the universe. We showed that, in both cases, the temperature remains frozen for a period 
at a plateau value TpL, which turns out to be much lower than its maximal value achieved 
during a pure reheating with the same initial (^energy density. 

As regards the ^^h? computation, we discriminated two basic types of x-production 
depending whether x's do or do not reach chemical equilibrium with plasma. In the latter 
case, two subcases were singled out: the type I and type II non-EP. The type I non-EP is 
activated for TpL ^ m-^^jlQ and AT,^ 7^ is required so as sizable Oj^/i^ is achieved. The 
type II non-EP is activated for TpL ~ m^/20. Finally, EP is applicable for TpL S> m^/2{). 

Next, we investigated the dependence of Q^h? on the ^^(Tj^g) variations, generated by 
varying the free parameters (p^j, m^, T^). We showed that mostly Q^h'^ increases with 
J7q(Tj^g) when TpL increases and it decreases as Qg(Tj^g) increases when TpL decreases, too. 
Also, Q^hP' decreases as {av) increases for EP and non-EPI for large (<Ju)'s, it decreases 
with {av) for non-EPII and large {av)^s and it remains constant for non-EPI and non-EPII 
for low {av)^s. Finally, in any case, ^x^"^ increases with and rUx- 

Comparing the results on O^/i^ with those in the QKS and the LRS, we concluded 
that in the present scenario, Q^h'^ does not exclusively increases with Q,q{Tj^g) (in contrast 
with the QKS) and it approaches its value in the LRS as Q,q{T^g) decreases. Finally, 
regions consistent with the present CDM bounds were constructed, using m^s and {av)^s 
commonly allowed in several particle models. In most cases, the required TpL is lower than 
about 40 GeV. As a consequence, simple, elegant and restrictive particle models such as 
the CMSSM [10] - which, due to the large predicted Jl^/i^, is tightly constrained in the SC 
[11] or almost excluded in the QKS [19, 34, 15] - can become perfectly viable in the KRS. 

Acknowledgments 

The author would like to thank K. Dimopoulos, G. Lazarides and A. Masiero for 
enlightening communications, I.N.R. Peddie for linguistic suggestions and the Greek State 
Scholarship Foundation (I. K. Y.) for financial support. 



References 



28 



References 

[1] D.N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003) [astro-ph/0302209]. 

[2] M. Tegmark et al, Phys. Rev. D 69, 103501 (2004) [astro-ph/0310723]; 
A.G. RiESS et al., Astrophys. J. 607, 665 (2004) [astro-ph/0402512] . 

[3] For a review from the viewpoint of particle physics, see 

A.B. Lahanas et al., Int. J. Mod. Phys. D 12, 1529 (2003) [hep-ph/0308251] . 

[4] For reviews, see K. Matchev, hep-ph/0402088; 
E.A. Baltz, astro-ph/0412170; 

G. Lazarides, hep-ph/0601016. 

[5] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983); 
J.R. Ellis et al., Nucl. Phys. B238, 453 (1984). 

[6] G. Servant and T.M.P. Tait, Nucl. Phys. B650, 391 (2003) [hep-ph/0206071]; 

H. C. Cheng et al, Phys. Rev. Lett. 89, 211301 (2002) [hep-ph/0207125]; 

K. Agashe and G. Servant, Phys. Rev. Lett. 93, 231805 (2004) [hep-ph/0403143]; 
J.A.R. Cembranos et al, Phys. Rev. Lett. 90, 241301 (2003) [hep-pli/0302041] . 

[7] E.W. KOLB AND M.S. Turner, The Early Universe, Redwood City, USA: Addison-Wesley (1990). 

[8] N. Okada AND O. Seto, Phys. Rev. D 70, 083531 (2004) [hep-ph/0407092] ; 
T. NiHEi, N. Okada and O. Seto, hep-ph/0409219. 

[9] M.S. Turner, Phys. Rev. D 33, 889 (1986); 

L. Covi et al., J. High Energy Phys. 06, 003 (2004) [hep-ph/0402240] ; 
J. Ellis et al., Phys. Lett. B 588, 7 (2004) [hep-ph/0312262] . 

[10] G.L. Kane et al., Phys. Rev. D 49, 6173 (1994) [hep-pli/9312272] . 

[11] J. Ellis et al., Phys. Lett. B 565, 176 (2003) [hep-ph/0303043] ; 

H. Baer and C. Balazs, J. Cosmol. Astropart. Phys. 05, 006 (2003) [hep-ph/0303114]; 
A.B. Lahanas and D.V. Nanopoulos, Phys. Lett. B 568, 55 (2003) [hep- ph/0303130]; 
U. Chattopadhyay et al., Phys. Rev. D 68, 035005 (2003) [hep-pli/0303201]. 

[12] U. Chattopadhyay and D.P. Roy, Phys. Rev. D 68, 033010 (2003) [hep-ph/0304108] . 

[13] T. Gherghetta, G.F. Giudice and J.D. Wells, JVucJ. Phys. B559, 27 (1999) [hep-ph/9904378] . 

[14] C. Pallis, Astropart. Phys. 21, 689 (2004) [hep-ph/0402033] . 

[15] C. Pallis, J. Cosmol. Astropart. Phys. 10, 015 (2005) [hep-ph/0503080] . 

[16] M. Kamionkowski and M.S. Turner, Phys. Rev. D 42, 3310 (1990). 

[17] J. McDonald, Phys. Rev. D 43, 1063 (1991); 

T. Nagano and M. Yamaguchi, Phys. Lett. B 438, 267 (1998) [hep-ph/9805204] . 

[18] G.F. Giudice, E.W. Kolb and A. Riotto, Phys. Rev. D 64, 023508 (2001) [hep-ph/0005123|. 

[19] P. Salati, Phys. Lett. B 571, 121 (2003) [astro-ph/0207396] . 

[20] R. Catena et al, Phys. Rev. D 70, 063519 (2004) [astro-ph/0403614] . 

[21] R. Allaiiverdi and M. Drees, Phys. Rev. Lett. 89, 091302 (2002) [hep-ph/0203118]. 

[22] R.J. Sciierrer and M.S. Turner, Phys. Rev. D 31, 681 (1985). 

[23] N. FORNENGO, A. Riotto and S. Scopel, Phys. Rev. D 67, 023514 (2003) [hep-ph/0208072] . 

[24] T. MOROI AND L. Randall, JVucJ. Phys. B570, 455 (2000) [hep-ph/9906527] . 

[25] M. FUJII AND K. Hamaguchi, Phys. Rev. D 66, 083501 (2002) [hep-ph/0205044] ; 
M. FUJII AND M. Ibe, Phys. Rev. D 69, 035006 (2004) [hep-ph/0308118]. 



29 



KINATION-DOMINATED REHEATING AND CDM ABUNDANCE 



[26] One more analysis of the LRS has been recently presented in ref. [27]. The authors used lower and 
higher (av) 's than those considered in ref. [14] and they applied the formalism to specific SUSY models. 
Nevertheless, our results in ref. [14] have been impressively verified - compare e.g. fig. 6 of the second 
paper (version of 1 May 2006) in ref. [27] with figs. 4-(ai) and 4-(a2) of ref. [14]. Moreover, in ref. [14] 
both thermal and non thermal contributions in the relevant Boltzmann equations have been taken 
into account contrary to what incorrectly mentioned in ref. [28]. 

[27] G. Gelmini and P. Gondolo, hep-ph/0602230: 

G. Gelmini, P. Gondolo, A. Soldatenko and C.E. Yaguna, hep-ph/0605016. 

[28] M. Drees, H. Iminniyaz and M. Kakizaki, Phys. Rev. D 73, 123502 (2006) [hep-ph/0603165]. 

[29] R.R. Caldwell et al, Phys. Rev. Lett. 80, 1582 (1998) [astro-ph/9708069]. 

[30] B. Spokoiny, Phys. Lett. B 315, 40 (1993) [gr-qc/9306008] ; 
M. Joyce, Phys. Rev. D 55, 1875 (1997) [hep-pli/9606223] ; 

P.G. Ferreira and M. Joyce, Phys. Rev. D 58, 023503 (1998) [astro-pli/9711102] . 
[31] C. Wetterich, JVucJ. Phys. B302, 668 (1988). 

[32] U. FRANgA and R. Rosenfeld, J. High Energy Phys. 10, 015 (2002) [astro-ph/0206194] . 
[33] C.L. Gardner, JVucJ. Phys. B707, 278 (2005) [astro-ph/0407604] . 

[34] S. Profumo and P. Ullio, J. Cosmol. Astropart. Phys. 11, 006 (2003) [hep-ph/0309220] . 

[35] A. Liddle and L.A Urena-Lopez, Phys. Rev. D 68, 043517 (2003) [astro-ph/0302054] . 

[36] B. Feng and M. Li, Phys. Lett. B 564, 169 (2003) [hep-ph/0212233]. 

[37] P.J. Peebles and A. Vilenkin, Phys. Rev. D 59, 063505 (1999) [astro-ph/9810509]; 
M. Yahiro et al, Phys. Rev. D 65, 063502 (2002) [astro-ph/0 106349]; 
K. Dimopoulos and J.W. Valle, Astropart. Phys. 18, 287 (2002) [astro-ph/0111417]; 
K. Dimopoulos, Phys. Rev. D 68, 123506 (2003) [astro-ph/0212264]. 

[38] E.J. COPELAND et al, Phys. Rev. D 64, 023509 (2001) [astro-ph/0006421]. 

[39] D.H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [hep-ph/0110002]. 

[40] E.W. KOLB, A. NOTARI AND A. RiOTTO, Phys. Rev. D 68, 123505 (2003) [hep-ph/0307241] . 

[41] Throughout this work we adopt the perturbative approach to the 4>-particle decay. However, if the 
initial amplitude of the 4>-oscillations is large enough [42], a stage of parametric resonance may take 
place, leading to explosive particle production known as preheating [43]. Applying this mechanism, a 
scenario similar to the one considered in this paper has been analyzed in ref. [44]. 

[42] L. KOFMAN, A. LiNDE AND A. Starobinsky, Phys. Rcv. Lett. 73, 3195 (1994) [hep-th/9405187]; 
L. KOFMAN, A. LiNDE AND A. Starobinsky, Phys. Rcv. D 56, 3258 (1997) [hep-ph/9704452] . 

[43] For a review, see B.A. Bassett, S. Tsujikawa and D. Wands, astro-ph/0507632. 

[44] M. Joyce and T. Prokopec, Phys. Rev. D 57, 6022 (1998) [hep-ph/9709320] . 

[45] D.J.H. Chung, E.W. Kolb and A. Riotto, Phys. Rev. D 60, 063504 (1999) [hep-ph/9809453] . 

[46] G. Belanger et al., Comput. Phys. Commun. 149, 103 (2002) [hep-ph/0112278] ; 
G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, hep-ph/0405253. 

[47] P. Gondolo et al., J. Cosmol. Astropart. Phys. 07, 008 (2004) [astro-ph/0406204] . 

[48] S. Hannestad, Phys. Rev. D 70, 043506 (2004) [astro-ph/0403291]. 

[49] R.H. Cyburt et al., Astropart. Phys. 23, 313 (2005) [astro-ph/0408033] . 

[50] R. Bean, S.H. Hansen and A. Melchiorri, Phys. Rev. D 64, 103508 (2001) [astro-ph/0104162]; 
JVucJ. PJiys. 110, (Proc. Suppl.) 167 (2002) [astro-ph/0201127]. 



References 



30 



[51] M.R. DE Garcia Maia, Phys. Rev. D 48, 647 (1993); 

M. Giovannini, Phys. Rev. D 60, 123511 (1999) [astro-pli/9903004] ; 

V. Sahni, M. Sami and T. Souradeep, Phys. Rev. D 65, 023518 (2002) [gr-qc/0105121]. 
[52] C.L. Bennett et al, Astropart. J. 464, LI (1996) [astro-ph/9601067]. 
[53] S. Davidson and S. Sarkar, J. High Energy Phys. 11, 012 (2000) [hep-pli/0009078] . 

[54] Assuming no entropy production, analytical expressions for this part of the rix evolution, have been 

recently presented in ref. [28]. 

[55] For a review, see C. MuNOZ, Int. J. Mod. Phys. A 19, 3093 (2004) [hep-ph/0309346l . 
[56] P. GONDOLO AND G. Gelmini, Nucl. Phys. B360, 145 (1991). 

[57] J. Ellis et al., Astropart. Phys. 13, 181 (2000) (E) ibid. 15, 413 (2001) [liep-ph/9905481]; 

M.E. Gomez, G. Lazarides and C. Pallis, Phys. Rev. D 61, 123512 (2000) [hep-ph/9907261]. 

[58] J. Edsjo and p. Gondolo, Phys. Rev. D 56, 1879 (1997) [hep-ph/9704361]. 

[59] C. BCEHM, A. Djouadi and M. Drees, Phys. Rev. D 62, 035012 (2000) [hep-ph/99 11496]; 

J. Ellis, K. Olive and Y. Santoso, Astropart. Phys. 18, 395 (2003) [liep-ph/0112113]; 
C. Pallis, Nucl. Phys. B678, 398 (2004) [hep-ph/0304047] . 

[60] A.B. Lahanas et al, Phys. Rev. D 62, 023515 (2000) [hep-ph/9909497] ; 
J. Ellis et al, Phys. Lett. B 510, 236 (2001) [hep-ph/0102098]; 

M.E. Gomez, G. Lazarides and C. Pallis, JVucJ. Phys. B638, 165 (2002) [hep-ph/0203131] 
{for an update, see G. Lazarides and C. Pallis, hep-ph/0406081). 



